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Newtonian and Post Newtonian Approaches 
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Regular and chaotic test particle motion in gravitational fields due to stellar bodies with quadrupo- 
lar and octupolar deformation are studied using Poincare surfaces of section. In first instance, we 
analyze the purely Newtonian case and we find that the octupolar term induces a distortion in the 
KAM curves corresponding to regular trajectories as well as an increase in chaoticity, even in the 
case corresponding to oblate deformation. Then we examine the effect of the first general relativistic 
corrections, provided by the post Newtonian approach. For typical values of the post Newtonian 
multipoles we find that the phase-space structure practically remains the same as in the classical 
case, whereas that for certain larger values of these multipoles the chaoticity vanishes. This impor- 
tant fact provides an interesting example of a situation where a non-integrable dynamical system 
. becomes integrable through the introduction of a large perturbation. 
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INTRODUCTION 



As is suggested by a wide variety of observational evidences, many astrophysical objects can be modeled as axially 
symmetric bodies with prolate or oblate deformation. For example, it is known that the Earth has non vanishing 
quadrupolar and octupolar moments, as a consequence of its oblate shape Q]. Also, many galaxies with a large disc 
component can be assumed as axisymmetric oblate bodies with a large quadrupolar moment and, in some cases, 
with a comparable octupolar deformation due to the remaining components (for example, the halo). Likewise, there 
are galaxy clusters with a cigar like shape and many dwarf galaxies that can be considered as nearly axisymmetric 
prolate deformed objects^. Although in the above cases the quadrupolar moment it is considered to be the major 
deviation from the spherical symmetry, there are situations (for example, some metallic clusters) where the octupolar 
[ deformation play a significant role0,Q. 

The motion of test particles around such stellar objects is a problem of wide physical interest. The case of attraction 
centers described by monopolar plus quadrupolar terms has been extensively studied by Gueron and Letelier, from a 
classical and relativistic standpoint, showing that the inclusion of external multipolar moments can induce chaos 0,0, 
(TlT | . In Newtonian gravity, as well as in general relativity, chaos can be found when the source has prolate deformation 
and the chaoticity grows by increasing the quadrupolar moment. On the other hand, it seems to be that the case 
corresponding to oblate deformation does not lead to chaotic motion, indeed for a very large quadrupolar deformation. 
However, inclusion of octupolar deformation can also induce chaos, as it was shown by Heiss, Nasmitdinov and 
RadufToj and Li [3] in the case of an harmonic oscillator. 

In the present paper we investigate the motion of test particles in axially symmetric gravitational potentials that 
are the sum of a monopole, a quadrupole and an octupolar term (only external multipole moments will be considered). 
At the first instance, in Sec. HH we shall perform the analysis in the context of Newtonian gravity by examining how 
the structure of the Poincare surfaces of section is determined by the octupolar moment. We find that, even in the 
case of oblate deformation, modest values of the octupolar moment induce chaotic motion. Then, in order to compare 
how these results are modified by considering the first general relativistic corrections, we shall introduce in Sec. IIIII 
the post Newtonian (1PN) approach, where the gravity is described by the classical Newtonian potential and another 
additional field, which is written in terms of the 1PN multipolar moments. Contrary to the usual statement that 
integrable problems in Newtonian theory become non integrable in general relativity, we find that, for certain values 
of the 1PN multipolar moments, the test particle motion is regular in situations where the Newtonian theory predicts 
chaos. 
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II. TRANSITION RING-SPINDLE TORUS AND CHAOS INDUCED BY OCTUPOLAR 

DEFORMATION 



Consider a test particle moving in the axially symmetric gravitational field generated by a stellar body with 
quadrupolar and octupolar deformation. In cylindrical coordinates (R,z,(p), its potential has the form 

a (3{2z 2 -R 2 ) j{2z 3 -3zR 2 ) 

~~VR 2 + z 2 ~ 2{R 2 + z 2 fl 2 ~ 2{R 2 + *2)7/2 • W 

Since we are interested on to describe the exterior test particle motion, we consider only external multipolar moments. 
In the above equation a is the monopole term that equals to Gm, where m is the total mass of the source and G is the 
gravitational constant. The quadrupolar term, denoted by /?, usually represents the major deviation from spherical 
symmetry In particular, if (3 > the source has prolate deformation and if (3 < we have the case corresponding to 
oblate deformation. The octupolar moment 7 describes the asymmetry of the source with respect to the equatorial 
plane, i.e. its "shape-of-pear" deformation. Both (3 and 7 are related to the source's density p(R, z) through the 
equations [9j 

/>OC />TT 

(3 = 2ir / r' 4 dr' / dd' sin 8' P 2 (cos 8') p(r' , 8'), (2) 
Jo Jo 

7 = 2tt / r' 5 dr' / d8' sin 8' P 3 (cos 8') p{r' ,8'), (3) 
Jo Jo 

where we have used spherical coordinates r = V 1 R 2 + z 2 , cos 9 = zjsj R? + z 2 and Pi denotes the Legendre polinomial 
of order I. 

The motion of a test particle in a gravitational field described by (Q]) , obeys the relations 

R = V R , (4) 

z 



V„ (5) 
dR ' 



where $ e / is the effective potential, given by 



*■/ = * + j^- W 



Here, £ = R 2 ip is the axial specific angular momentum that is conserved as a consequence of the axial symmetry. The 
second integral of motion is the total specific energy, given by 

E=(V 2 + V 2 )/2 + <S> e} . (9) 

According to eqs. (@|-((9|), the motion is restricted to a three dimensional phase space (R, z, Vr). This fact enable us 
to introduce the Poincare surfaces of section method, in order to investigate the trajectories of test particles. 

In Fig. [T] we plot a typical z — surface of section corresponding to test particle motion in the presence of a 
gravitational field due to a prolate deformed source whom octupolar moment vanishes (from here on, we choose 
a = 1, without loss of generality). We note a central and lateral regular regions composed by ring torus (ring KAM 
curves). They are enclosed by a chaotic region containing two small resonant islands near its top and bottom edges. 
In Fig. [5] we switch on the octupolar moment (7 = 0.02), maintaining the same prolate deformation. The resulting 
surface of section presents a more prominent chaotic region, since in this case the outer zones of resonant islands have 
overlapped. The regular regions now contains "spindle" torus. They can be viewed clearly in the central region and 
scarcely insinuated in the lateral zone. In Fig. [31 as a consequence of increase the octupolar moment to 7 = 0.04, 
the chaotic region is more prominent (the lateral regular zone has disappeared), as well as the central spindle KAM 
curves. Finally, Figs. [4] and [5] show the effect caused by the progressive rise in the octupolar deformation, starting 
from a regular prolate situation. 

The fact that by switching on the octupolar moment increases the chaoticity and leads to apparition of spindle 
torus, can be viewed even in the case corresponding to oblate deformation, which commonly presents regular motion 
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(in fact, there is a wide numerical evidence that particles moving around a monopole plus an oblate quadrupole are 
not chaotic). Fig. [5] shows the transition from regularity to chaos by increasing 7. With E = —0.32, t = 1.1 and 
(3 = —0.2, we start from 7 = (Fig. [6^ ) which corresponds to a regular motion. The case 7 = 0.02 is also regular but 
the KAM curves has been distorted to spindles. In Fig. [5J; (7 — 0.04) the distortion is more prominent and finally, 
when 7 has increased up to 0.06 we note the apparition of a chaotic region enclosing the spindles. 



III. NEWTONIAN VS POST-NEWTONIAN MULTIPOLAR MOMENTS 

The post newtonian approximation (commonly called as 1PN approach) gives the first general relativistic corrections 
(up to order v 2 /c 2 , where c is the speed of light) to the motion equations, when we deal with particles moving non 
relativistically (v <C c) in the presence of strong gravitational fields. We are interested in the situation when a 
test particle moves around a static axially symmetric source. As it was shown by Rezania and Sobouti (l3| . in the 
stationary case we have an integral of motion E = v 2 /2 + $ + (2<I> 2 + ty)/c 2 that can be considered as a generalization 
of the classical specific energy. Here, 5> and \1/ obeys the Poisson equations 

V 2 $ = ^t 00 , (10) 
2 2 

V 2 * = 47rG[T 00 + T% (11) 

02 2 
where T 00 and T 00 denotes the 0-th and 2-th terms in v/c in the 00-components of the energy momentum tensor. T n 
represents the trace of T y up to order 2 in v/c. Relations (jTDJ) and (fTTj) suggest that both $ and can be expanded 



in external multipolar fields. In particular, since T 00 represents the source's density [l2l [l3|. the expansion for $ is 
the same as in equation (|TJ) - 
In the case of we can write 

, T ,_ & P(2z 2 -R 2 ) 7 (2z 3 -3zfl 2 ) 

Vi? 2 + z 2 2(i? 2 + z 2 ) 5 /2 2 (i? 2 + z 2 )V2 ' 1 > 

where a, (3, and 7 play the role of the post Newtonian monopole, quadrupole and octupole moments, respectively. In 
a similar way as in ([l])-©, they are given by 

/•oo pi 2 2 

a = 2it J r' 2 dr' J dx'[T m + T% (13) 

/.oo pi 2 2 

= ^J r' 4 dr' J dx'P 2 (x')[T 00 + T% (14) 

/•oo pi 2 2 

7 = 2tt / r' 5 dr' / dx' P 3 (x')[T 00 + T% (15) 



where x' = cos 8' = z'/V R' 2 + z' 2 . 

The motion equations for a test particle can be derived from the Hamiltonian 



H= 1 -(V 2 + V z 2 ) + ^ + <S> + ±(2<S> 2 + *), (16) 



(note that I — R tp is an integral of motion) and can be written as 



where we have defined the post Newtonian effective potential U e f as 



V, - (13) 



t/ ./ = ' |, + ^( 24,2 + *) + 2l5' (19) 
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From eqs. (fl7|) . (fT8|) and (|19p . along with ([T]) and (fT2|) , we can study the trajectories of test particles in a similar way 
than in the Newtonian case, as was performed in section HT1 

For example, Fig. [8] show a, z — surface of section corresponding to I = 0.9, E = —0.4, a = 1, = 0.3, 7 = 0.02 
(the same parameters as in Fig. [3]), a — 10 5 , (3 = —2 x 10 4 and 7 = 10 4 (according with the system of units employed, 
we have to take c = 10 4 ). It has similar features as in Fig. [3l but in this case there is not lateral resonant islands. If 
one increase adequately the post Newtonian quadrupole and octupole, an interesting fact occurs: the chaotic zones 
disappear and the resulting surface of section reveals completely regular motion. Indeed, as it is shown in Fig. [9J 
the regular regions are made by rings instead of spindle torus, for these particular values of multipole parameters. In 
Figure [7] we show the contours of the effective potential corresponding to these two cases. The former exhibits z = 
asymmetry while the second is practically symmetric with respect to the equatorial plane. 



IV. CONCLUDING REMARKS 



Octupolar deformation in astrophysical objects can introduce significant modifications to the phase-space structure 
corresponding to test particles moving around prolate or oblate centers of attraction. Apart from an increasing in 
the chaoticity, the apparition of spindle torus in regular regions is a remarkable effect caused by the asymmetry of 
the source with respect to its equatorial plane. A larger equatorial asymmetry involves more distorted spindle KAM 
curves and more prominent stochastic regions in the phase-space. This fact carries dramatic consequences in the case 
of oblate deformed sources, which usually are associated with regular motion. Here, chaos emerges once we switch on 
the octupole moment. 

When the first general relativistic corrections are taken into account, through the post Newtonian approach, the 
structure in the phase space remains practically the same as in Newtonian gravity, for typical values of 1PN multipolar 
moments. However, there are certain large values for which the test particle motion becomes regular. This fact intro- 
duces an interesting example where a non-integrable dynamical system becomes integrable through the introduction 
of a large perturbation. 
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FIG. 2: Surface of section for the same initial conditions of the previous figure. We maintain the values i — 0.9, E = —0.4, 
a = 1 and (3 = 0.3, but now 7 = 0.02. 
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FIG. 3: In this case, setting the same values of the two previous figures but 7 = 0.04, we have a prominent chaotic zone 
enclosing three "spindle" KAM curves. 
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FIG. 4: Surfaces of section for I = 0.9, E = -0.4, a = 1, (3 = 0.2 and (a) 7 = 0; (b) 7 = 0.02; (c) 7 = 0.04; (d) 7 = 0.06. In 
each case, they are generated by three orbits with initial conditions (i) z = 0, R = 0.91, Vr = 0; (ii) z = 0, R = 0.78, Vr = 
and (iii) z = 0, R= 1.16, Vr = 0.18. 
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FIG. 5: An orbit in the meridional plane with initial conditions z = 0, R = 0.78, Vr = and the same parameters considered 
above. Again we have the cases (a) 7 = 0; (b) 7 = 0.02; (c) 7 = 0.04; (d) 7 = 0.06 




FIG. 6: Surfaces of section for £ = 1.1, E = -0.32, a = 1, /3 = -0.2 and (a) 7 = 0; (b) 7 = 0.02; (c) 7 = 0.04; (d) 7 = 0.06. In 
each case, they are generated by three orbits with initial conditions (i) z = 0, R = 0.91, Vr = 0; (ii) z = 0, R = 0.78, Vr — 
and (iii) z = 0, R = 1.16, Vr = 0.18. 
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FIG. 9: Now we keep the same values of E, I, a, f3, 7, a as well as the same initial conditions as in previous figure, but 
P — — 4 x 10 7 and 7 = -2 x 10 6 . With these values of post Newtonian multipolar moments we have a surface of section 
corresponding to regular motion 



